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NNLO r+r production cross section at threshold 
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The threshold behaviour of the cross section a{e'^e~ — > r^r~) is analysed, taking into account the known 
higher-order corrections. At present, this observable can be determined to next-to-next-to-leading order (NNLO) 
in a combined expansion in powers of Os and fermion velocities. 



1. Introduction 



The Tau-Charm Factory, a high-luminosity 



(~ 10' 
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collider with a centre- 



of-mass energy near the r+r" production thresh- 
old, has been proposed as a powerful tool to per- 
form high-precision studies of the r lepton, charm 
hadrons and the charmonium system In re- 
cent years, this energy region has been only par- 
tially explored by the Chinese BEBC machine 
(~ 10^^ cm~^ The possibility to operate the 

Cornell CESR collider around the t+t^ thresh- 
old has revived again the interest on Tau-Charm 
Factory physics. 

A precise understanding of the e+e^ t^t^ 
production cross section near threshold is clearly 
required. The accurate experimental analysis of 
this observable could allow to improve the present 
measurement Q of the r lepton mass. The cross 
section a{e~^e~ —> t'^t^) has already been anal- 
ysed to 0{a'^) in refs. [Q, including a resumma- 
tion of the leading Coulomb corrections. 

The recent development of non-relativistic ef- 
fective field theories of QED (NRQED) and QCD 
(NRQCD) m has allowed an extensive investiga- 
tion of the threshold production of heavy flavours 
at e+e~ colliders. The threshold bb and tt 
|6| production cross sections have been computed 
to the next-to-next-to-leading order (NNLO) in 
a combined expansion in powers of and the 
fermion velocities. Making appropriate changes, 
those calculations can be easily applied to the 
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study of T+T production 



One can then 



achieve a theoretical precision better than 0.1%. 

2. Perturbative calculation to 0{a'^) 

At lowest order in QED, the r leptons are pro- 
duced by one-photon exchange in the s-channel, 
and the total cross section formula reads 



(1) 



where v ~ -y/l — 4to^ / s is the velocity of the final 
r leptons in the center-of-mass frame of the e+ e~ 
pair which makes as vanish when v —^ 0. 

Electromagnetic corrections of 0{a) arise from 
the interference between the tree level result and 
1-loop amplitudes. A factor a/v emerges in the 
1-loop final state interaction between the tau lep- 
tons, making the cross section at threshold finite. 
Furry's theorem guarantees that contributions to 



a(e^ 



') coming from initial, intermedi- 



ate and final state corrections completely factor- 
ize at 0{a^), including real photon emission. 

Some undesirable features appear at 0{a'^): 
The two-loop T+T^7 vertex develops an c? jv^ 
term which makes the cross section ill-defined 
when w ^ 0, and multiple photon production of 
tau leptons by box-type diagrams and the non- 
zero interference of initial and final state radition 
spoil exact factorization. However, as it has been 
recently shown [H, the squared amplitude of the 
e^e~ — > T+T" box diagram is proportional to 
OL^v^ , and so represents a N'^LO correction in the 
combined expansion in powers of a and u, far 
beyond the scope of this analysis. In addition. 



1 



contributions to the total cross section from di- 
agrams with real photons emitted from the pro- 
duced taus can be shown to begin at N'^LO, and 
factorization remains at NNLO. The total cross 
section can thus be written as an integration over 
the product of separate pieces including initial, 
intermediate and final state corrections: 



F{s,w) 



1 



1 -I- e2no.„(u)) 



a{w)dw. (2) 



The radiation function F{s,w) describes ini- 
tial state radiation, including virtual corrections. 
The integration accounts for the effective energy 
loss due to photon emission from the e"*" e~ pair, 
and it includes the largest corrections coming 
from the emission of an arbitrary number of ini- 
tial photons, which can sizeably suppress the total 
cross section, a collects only final-state interac- 
tions between the tau leptons, and it is usually 
written in terms of the tau spectral density R^, 



T+T-)=R^{s)a, 



pt •) 



(3) 



with 



4-n-a 
3s 



The threshold behaviour of the 



total cross section will be ruled by the expansion 
of at low velocities. 

3. Non-Relativistic Corrections: NRQED 

A NNLO calculation of the cross section in 
the kinematic region where v ^ a has to ac- 
count for all terms proportional to v{a/v)" x 
[l;a;v;a'^;av;v'^] with n = 1,2,... The leading 
divergences (i.e. (^) , n > 1) are resummed in 
the well-known Sommerfeld factor iflOl 



I*c,b(0)|^ 



1 — exp(— QfTr/t;) 



(4) 



multiplying the Born cross section (|l|). A sys- 
tematic way to calculate higher-order corrections 
in this regime requires the use of a simplified 
theory which keeps the relevant physics at the 
scale Mv ~ Ma, characteristic of the Coulomb 
interaction. NRQED jj] was designed precisely 
for this purpose. It is an effective field theory 
of QED at low energies, applicable to fermions 
in non-relativistic regimes, i.e. with typical mo- 
menta p/M ~ V <C 1. Interactions contained in 



the NRQED Lagrangian have a definite velocity 
counting but propagators and loop integrations 
can also generate powers of v. With appropriate 
counting rules at hand, one can prove that all in- 
teractions between the non-relativistic pair t+t~ 
can be described up to NNLO in terms of time- 
independent potentials , derived from the low- 
energy Lagrangian. Therefore, the low-energy ex- 
pression of the r spectral density is related with 
the non-relativistic Green's functions [[12|: 



R": 



6 TT 



Im(Ci G(i?) ~|^G,(i?)) (5) 



with Ci a short distance coefficient to be deter- 
mined by matching full and effective theory re- 



sults and E — 



the non-relativistic energy. 



The details of this derivation can be found in the 
Appendix B of 0. The Green's function G obeys 
the Schrodinger equation corresponding to a two- 
body system interacting through potentials de- 
rived from £nrqed at NNLO: 



+ Vc{r) + V^^{r) + VUr)-E 



M AAP 

X G{r,r\E) = S^^^r - r') 



(6) 



The term 



4M' 



is the first relativistic correction 



to the kinetic energy. Vc stands for the Coulomb 
potential with 0(0^) corrections. At NNLO, the 
heavy leptons are only produced in triplet S-wave 
states, so we just need to consider the corre- 
sponding projection of the Breit-Fermi potential 
Vbf. Finally, Vao is a NNLO piece derived from 
a contact term in £nrqed, which reproduces the 
QED tree level s-channel diagram for the process 

A solution of eq. must rely on numerical 
or perturbative techniques. In the QED case, 
a significant difference between both approaches 
is not expected, being a such a small parame- 
ter. Consequently we follow the perturbative ap- 
proach, using recent results where NLO and 
NNLO corrections to the Green's function are 
calculated analytically, via Rayleigh-Schrodinger 
time- independent perturbation theory around the 
known LO Coulomb Green's function Gc- We re- 
fer the reader to Appendix C of for complete 
expressions of the Green's function corrections. 
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4. Vacuum Polarization 

For a complete NNLO description of (T(e+e~ 
T~^T~), two-loop corrections to the photon prop- 
agator should be included. The light lepton con- 
tributions to the vacuum polarization are the 
standard 1- and 2-loop perturbative expressions. 
For the r contribution in the threshold vicinity 

> 4M^, resummation of singular terms in the 
limit u ^ is again mandatory. Under the as- 
sumption a ^ V, we need to know NLO contribu- 
tions to Ylriq^), performing the direct matching 
for both real and imaginary parts. 

In the hadronic sector, we can relate the 
hadronic vacuum polarization with the total cross 
section a{e^e^ ^ ^ had). Below 1 GeV, the 
electromagnetic production of hadrons is domi- 
nated by the p resonance and its decay to two 
charged pions. The photon mediated tt+tt" pro- 
duction cross section is driven by the pion elec- 
tromagnetic form factor F(s). An analytic ex- 
pression for F{s) was obtained in Ref. using 
Resonance Chiral Theory and the restrictions im- 
posed by analyticity and unitarity. The obtained 
F(s) provides an excellent description of exper- 
imental data up to energies of the order Sp ^ 
1 GeV^. 

For the integration region above Sp, we use 
Imlltad as calculated from pQCD. Our simple es- 
timate has been proved to deviate by less than 
5% for the running of a at the scale y/s = Mz- 
Considering that lii^^^ modifies cr{e^e~ r+r^) 
near threshold by roughly 1%, our result has a 
global uncertainty smaller than 0.1% for the to- 
tal cross section. 

5. Bound states effect 

Green's functions develop energy poles below 
threshold corresponding to spin triplet (n'^S'i) 
electromagnetic t+t" bound states. The small 
width of these bound levels, dominated by their 
e+e" decay rate 

rrira^ 6.1 • IQ-^ eV 
Fee = -r^ « 5 , (7) 

make these states very difficult to be resolved ex- 
perimentally. For the same reason, the bound 
states cannot affect the shape of the cross section 
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Figure 1. The spectral density Rr at low veloci- 
ties in both QED and NRQED. 

at threshold, contrary to the case of heavy quark 
threshold production, where bound states play a 
crucial role. 

6. Numerical analysis for a{e^e~ r+r ) 

The need for performing resumma- 
tions of the leading non-relativistic terms 
(a/v)" [v, va, ^ . . .] is evidenced in Fig. ||. The 
QED spectral density vanishes as w ^ 0, due to 
the phase space velocity in formula (|l|), which 
is cancelled by the first term appearing in 
the O{o) correction, making the cross section 
at threshold finite. More singular terms near 
threshold, w^'^, . . . arising in higher-order correc- 
tions completely spoil the expected good conver- 
gence of the QED perturbative series at low v. 
This is no longer the case for the effective theory 
perturbative series, whose convergence improves 
as we approach the threshold point, as shown in 
Fig. 1^, and higher-order corrections reduce the 
perturbative uncertainty. In the whole energy 
range displayed in Fig. |^, the differences between 
the NNLO, NLO and LO results are below 0.8%, 
which indicates that the LO result, i.e. the Som- 
merfeld factor, contains the relevant physics to 
describe the threshold region. Adding the inter- 
mediate and initial state corrections we have a 
complete description of the total cross section of 
r+T~ production, as shown in Fig. ^. Coulomb 
interaction between the produced r's, becomes 
essential within few MeV's above the threshold. 
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Figure 2. Relative sizes of corrections to Rt{s) as 
calculated in NRQED. 

and its effects have to be taken into account to 
all orders. Initial state radiation effectively re- 
duces the available center-of-mass energy for r 
production, lowering in this way the total cross 
section. 

We should emphasize that NNLO corrections 
do not modify the predicted behaviour of the 
LO and NLO cross section, but are essential to 
guarantee that the truncated perturbative series 
at NLO gets small corrections from higher-order 
terms. Hence, we have shown that the theoretical 
uncertainty of our analysis of cr(e+e^ t^t^) 
is lower than 0.1%. Nevertheless, the statistical 
uncertainty of the most recent experiments is still 
much larger than the theoretical one due to low 
statistics, and we should wait for future machines 
to improve it. 
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